Can pulling cause right- to left-handed structural transitions in negatively supercoiled 

DNA double-helix? 



Zhou Haijur£il 

Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, China 
State Key Lab. of Scientific and Engineering Computing, Beijing 100080, China 



CU- 
(N 



Ou-Yang Zhong-can 

Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, China 

(February 1, 2008) 

The folding angle distribution of stretched and negatively supercoiled DNA double-helix is investi- 
gated based on a previously proposed model of double-stranded biopolymers (H. Zhou et al, Phys. 
Rev. Lett. 82, 4560 (1999)). It is shown that pulling can transit a negatively supercoiled DNA 
double-helix from the right-handed B-form to a left-handed configuration which resembles DNA 
Z-form in some important respects. The energetics of this possible transition is calculated and the 
comparison with recent experimental observations is qualitatively discussed. 
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Because of its vital biological significance, the elasticity 
of DNA has invoked considerable interest during recent 
years, and it is now known experimentally that radical 
transitions in DNA internal structure can be induced by 
the action of mechanical forces and/or torques. For ex- 
ample, pulling a DNA chain with a force of 70 piconew- 
tons (pN) will convert DNA standard B-form conforma- 
tion into a over-stretched S-forrroo; and at the joint ac- 
tion of a positive torque and a force about 3 pN,_a DNA 
will take on a novel P-form with exposed basesa. Here 
we suggest a further possibility for_this transition sce- 
nario and show that under-twistedo (negatively super- 
coiled) DNA can take on a left-handed configuration un 
der the action of a moderate stretching force. Such a left 
handed conformation is found to resemble DNA Z-for 
in some important respects. The energetics of this possi- 
ble transition is calculated and-,a qualitative comparison 
with very recent experimentsElQ is also performed. 

The extension vs supercoihng relatione! for under- 
twisted DNA is studied based on a model proposed 
earlieiQ (see Fig. |j]a). For small pulling forces (< 0.3 
pN), a supercoiled DNA can shake off its twisting stress 
by writhing its central axis and forming plectonemic 
structures^. However, this leads to shortening of DNA 
end-to-end distance and hence becomes very unfavor- 
able as the force increases. At this stage, the torsional 
stress caused by supercoiling begins to unwind the B- 
form double-helix and triggers the transition of DNA in- 
ternal structure, where a continuously increasing portion 
of DNA takes on some certain new configuration as super- 
coiling increases, while its total extension keeps almost 
invariant. 

Information about the new configuration can be re- 
vealed by the folding angle! ip distributions P(y). This 
distribution is calculated by 



P(tp) = / <& 2 (t,^)dt, 



(1) 



where t is the tangent vector of DNA's central axis and 



$(t, ip) is the (normalized) ground-state eigenfunction of 
the Green equation Eq. (9) in Ref. 8. The folding 
angle distribution (Fig. 0b) has the following aspects: 
When the torsional stress is small (with the supercoiling 
degree \a\ < 0.025), the distribution has only one nar- 
row and steep peak at cp ~ +57.0°, indicating that DNA 
is completely in B-form. With the increase of torsional 
stress, however, another peak appears at ip ~ —48.6° 
and the total probability for the folding angle to be 
negatively-valued increases gradually with supercoiling. 
Since negative, folding angles correspond to left-handed 
configurationsa, we can conclude that, with the increas- 
ing of supercoiling, left-handed DNA conformation is nu- 
cleated and it then elongates along the DNA chain as 
B-DNA disappears gradually. The whole chain becomes 
completely left-handed at a ~ —1.85. 

It is worth noticing that, (i) as the supercoiling degree 
changes, the positions of the two peaks of the folding 
angle distribution remain almost fixed and, (ii) between 
these two peaks, there exists an extended region of fold- 
ing angle from to 7r/6 which always has only extremely 
small probability of occurrence. Thus, a negatively su- 
percoiled DNA can have two possible stable configura- 
tions, a right-handed B-form and a left-handed configu- 
ration with an average folding angle ~ —48.6°. A tran- 
sition between these two structures for a DNA segment 
will generally lead to an abrupt and finite variation in 
the folding angle. 

The sum of the average base-stacking energy and tor- 
sional energy caused by external torque (see caption of 
Fig. ||a) as a function of torque T is shown in Fig. |^a 
and the relation between a and T in Fig. From these 
figures we can infer that, (i) for negative torque less than 
the critical value T c ~ —3.8 k#T, DNA can only stay in 
B-form state; (ii) near r c DNA can either be right- or be 
left-handed and, as negative supercoiling increases (see 
Fig. |^b) more and more DNA segments will stay in the 
left-handed form, which is much lower in energy (~ —2.0 
k^T per base pair (bp)) but stable only when torque 
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reaches r c ; (iii) for negative torque greater than r c DNA 
is completely left-handed. B-form DNA at r c has energy 
about 0.0 k^T per bp, indicating that the work done 
by external torque just cancels the base-stacking energy. 
Therefore, it might not be enough to further break the 
hydrogen bonds between DNA complementary bases and 
cause denaturationoQ. Nevertheless, since the transition 
from right- to left-handed structure requires radical rear- 
rangement of DNA base pairs, the possibility of transient 
denaturation in DNA double-helix can not be ruled out. 
This is a subtle question, and maybe transient denatu- 
rations can occur in the weaker AT-rich regions, or even 
be induced and then captured by the addecLhomologous 
single-stranded DNA probes in the solutions. 

For the left-handed state revealed by Figs. [l]b and || 
we have obtained that, at T = —4.0 k^T (where DNA is 
completely left-handed) the average rise per bp is about 
3.83 A, and the pitch per turn of helix is 46.76 A, with 
the number of bps per turn of helix being 12.19. Notice 
these characteristic quantities are very similar with those 
of DNA left-handecLZ-form, which are 3.8 A, 45.6 A, and 
12 bps, respectively^. We suspect that the identified left- 
handed configuration belongs to DNA Z-form. Recently, 
Leger et al. also pointed out that Z-form structure should 
be included to qualitatively interpret their experimental 
resuhu. 



a ' Electronic mail: zhouhj@itp.ac.cn 
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FIG. 1. a, DNA extension 
percoiling degree a is denned 
and Ik are, respectively, the 1 
torsionally constrained DNAl 
are from Ref. 4. b, folding 
pulled with a force of 1.3 pN. 
fined in the range from — n/2 
values corresponding to left ( 



vs supercoiling degree. The su- 
as o~ = (Ik — lko)/lko, where lko 
king number of a relaxed and a 
Experimental data (symbols) 
angle distributions for a DNA 
Here, the folding angle y> is de- 
to 7r/2, with negative (positwe) 
right)-handed configurations!!!. 
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FIG. 2. a, The sum of average base-stacking and torsional 
energy per base pair at force 1.3 pN. For highly extended 
DNA only these two interactions are sensitivs-iwith torque. 
The torque-related torsional energy density isa — T sin ip/R, 
where R is DNA molecular radius, b, The relation between 
DNA supercoiling degree and external torque at force 1.3 pN. 



